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Abstract 

The temperature dependence of the pulse length was measured for an 18-fold seg- 
mented n-type germanium detector in the temperature range of 77 - 120 K. The 
interactions of 122 keV photons originating from a ^^^Eu source were selected and 
pulses as observed on the core and segment electrodes were studied. In both cases, 
the temperature dependence can be well described by a Boltzmann-like ansatz. 
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1 Introduction 

High purity germanium detectors, HPGeDs, are widely used in spectroscopy, 
for example in the detection of low levels of radioactivity. They are also used 
in a wide variety of applications in particle and nuclear physics [Tp] . In par- 
ticular, they are used for gamma ray tracking in arrays like AGATA [3] and 
GRETA |3], and in the search for neutrinoless double beta decay in experi- 
ments hke GERDA [1] and MAJORANA 0. 

For some applications, the shapes of the electric pulses collected on the elec- 
trodes are of interest. They are used in so called pulse shape analyses which 
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are usually used to obtain information about event topologies. This can either 
help in gamma ray tracking [6] or in identifying background events \7W\- 

A fundamental parameter of a pulse is its length. It can be used to measure the 
distance of the energy deposit from the electrodes. Its absolute value, however, 
depends also on the impurity density of the crystal and the relative location of 
the energy deposit with respect to the crystal axes. Other important factors 
are the mobilities, fi, of the charge carriers, electrons and holes, which are 
temperature, T, dependent. 

The mobilities of the charge carriers have been measured in the past. However, 
the results, obtained for different samples, do not completely agree [9J . In prin- 
ciple, the mobilities can be determined by comparing measured to simulated 
pulses. In practice, this is difficult, because the impurity density distribution 
of a detector crystal is also not precisely known. 

The T dependence of the mobilities, can be directly related to the T depen- 
dence of the pulse lengths. Photon interactions from the 122 keV line of a 
^^^Eu source were used to select charge depositions close to the outer surface 
of an n-type cylindrical true-coaxial germanium detector. Holes created close 
to the surface quickly drift to the mantle electrode and the electrons drifting 
to the core determine the total length of the pulse. The holes influence only 
the first part of the pulse. Pulses were studied separately for interactions close 
to the crystallographic axes (110) and (100). 



2 Experimental Setup 

The detector used was a cylindrical true-coaxial 18-fold (60, 3z) segmented 
n-type detector with a height of 70 mm and a diameter of 75 mm. It was 
manufactured by Canberra France and is similar to the detector previously 
characterized in detail [10]. The inner bore had a diameter of 10 mm. The 
density of electrically active impurities was given by the manufacturer as 0.35- 
10^°/cm^ at the top and 0.55 ■ 10^°/cm^ at the bottom of the detector. The 
change in impurities is assumed to be linear with height, z. The operational 
voltage was 2000 V. 

The detector was mounted either inside a vacuum cryostat [lOj, Kl, or sub- 
merged in a liquid nitrogen volume [TT], GIL In both collimated 40 kBq 
^^^Eu source with a la beam spot diameter of about 5 mm was used. The co- 
ordinate system used was cylindrical with the origin at the geometric center 
of the detector and = along its crystallographic (110) axis. 

The detector T inside GII was constantly 77.4 K while it was increasing with 
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time inside Kl, where it was cooled via a cooling finger submerged in a liquid 
nitrogen reservoir. When the reservoir was not refilled, the detector slowly 
warmed up. The temperature, T^^n, was monitored with a thermal resistor, 
PTIOO, mounted as close as possible to the detector inside Kl. The tempera- 
ture, T, of the detector was expected to be higher by AT = T — Tmon between 
4- lOK. 

The following datasets were analysed: 

• DSl: detector operated in GII, submerged in LN2, T = 77.4 K; 

• DS2: detector operated in Kl, vacuum, 95 K < Tmon < 100 K; 

• DS3: detector operated in Kl, vacuum, 100 K < Tmon < 120 K. 

The measurements in Kl were started after the initial cool-down of the setup. 
They were performed while the detector was warming up. For DS2, Tmon was 
monitored during breaks in the data taking. The value of Tmon is interpolated 
for intermediate times. For DS3, the time dependence of Tmon was taken from 
test warm-ups where Tmon was measured in one minute intervals. Figure [1] 
shows the temperature development for one measurement of DS2 and two test 
warm-ups for DS3, one at the beginning and one at the end of the DS3 data 
taking. 

The absolute value of Tmon for DS3 depended on the initial Tmon- The de- 
velopment of Tmon with time was however stable. The Tmon could thus be 
extrapolated from the initial Tmon for each measurement cycle. 




Fig. 1. Development of Tmon for one of the measurements of DS2 and two test 
warm-ups for DS3. 

Data were grouped into 15-minute intervals. For the analysis, it was assumed 
that T and Tmon were constant during these intervals. At the beginning of the 
measurements in DS3 the steep rise in Tmon was excluded. 
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Measurements were taken separately around the (110) axis at = 0° and 
the (100) axis at = 45°. The measurements for DS3 were 5° off axis. The 
restrictions of the Gil setup did not allow a measurement on the (100) axis. 
Additional data were taken for cross-checks. Table [T] provides an overview. 
The two axes were located in different segments of the detector. None of the 
axes was close to a segment boundary. A schematic of the setup is shown in 
Fig.H 
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Table 1 

Datasets Al and A2, entering the analysis, and CI and C2, used for cross checks 



^90° 

<110> 




Fig. 2. Schematic of the measurement setup. Long- and short-dashed lines indicate 
the segment boundaries and crystaHographic axes. 

The core and all 18 segments, 19 channels, were connected to charge sensi- 
tive preamplifiers, PSC 823, produced by Canberra-France. The preamplifiers 
and the cables limit the bandwidth of the system to about 10 MHz. The ener- 
gies and pulse forms were recorded for all 19 channels using a "DGF Pixie-4" 
data acquisition system, DAQ, manufactured by X-Ray Instrumentation As- 
sociates, XIA. The core pulse was always used to trigger. The DAQ has 14-bit 
ADCs and a sampling frequency of 75 MHz |12j . 
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3 Pulse shape simulation 



Two simulated pulses, one for the core and one for the segments, were used 
as references to determine the relative lengths and amplitudes of measured 
pulses. They were generated using the pulse shape simulation package [13] 
developed for true-coaxial detectors. 

The pulses were simulated for a point-like energy deposit of 122 keV at r = 
32.5 mm, = 0° and z = mm. This point is located on the (110) axis of 
the detector, 5 mm deep inside the detector. The depth was chosen as it is the 
average penetration depth of 122 keV photons. 

The parameters of the simulation were: 

• impurity level of 0.45 ■ 10^*^ /cm^, corresponding to the specification ai z = 
mm; 

• electron mobility [9j, corresponding to T =77 K; 

• grid for numerical field calculations of 32 (r) x 180(0) x 70 (z); 

• time, t, development in steps of 1 ns, corresponding to spatial steps of about 
100 fxm which is the granularity needed to reflect the development of charge 
trajectories; 

• bandwidth of 10 MHz and amplifier decay time of 50 /is; 

• no noise. 

The simulated charge pulses, Cg/^, for core, c, and segments, s, are shown in 
Fig. [3l The amplitude is normalized to ±1, respectively. The pulses start at 
t=100ns and have a total rise time, tr, of almost 400 ns. The times, at which 
the pulses reach 10 % and 90 % of their amplitudes, are indicated. The time 
interval inbetween, the so called 10-90 rise-times, ^^,5:10-90^ differ for core and 
segments. The t^,s:io-90 -^qj-q y^^ed in all measurements to evaluate the length 
of a pulse, because the beginning and the end of a pulse are very susceptible 
to noise. 

The usage of a reference pulse, simulated as one interaction 5 mm inside the 
detector, instead of an average over pulses simulated according to beamspot 
and penetration depth, does not change significantly [13]. The effect on 

j^s:io-90 ^jjj subject of a separate investigation. In general, the hole drift is 
less well understood than the electron drift. This is reflected in a less accurate 
prediction of the detailed pulse shape for the segments. 

The simulated core and segment pulses along (110) at 77.4 K have = 
290 ns and tf^-^° = 239 ns. 
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Fig. 3. Simulated reference pulses for a) core and b) segments. The pulses are nor- 
malized to ±1. The times at which 10 % and 90 % of the amplitude are reached are 
indicated by dashed lines. 



4 Model expectations 



The T dependence of the rise time is assumed to be due to the T dependence of 
the mobiUty of the charge carriers. The velocity, Ve/h{x) of the charge carrying 
electrons, e, and holes, h, at each point x = (r, 0, z) is given by 



'Je/h 



[x) 



fie/hE{x) 



[1 + 



■M2 



Eo 



(1) 



where fig/h is the charge carrier mobility tensor, E[x) is the electric field and 
Eq and /9 are fit parameters obtained in measurements [9] along the crystal- 
lographic axes. The mobility tensor depends on the relative position of x to 
the crystallographic axes. Along the axes, fie/h becomes a number. A more 
detailed discussion relevant for the detector considered here can be found 
elsewhere [T3] . 

As the detector is true-coaxial, E is radial with the absolute value, E{r). Thus, 
the velocity, Ve/h{r), becomes radial along the crystallographic axes. Due to 
the extremely low impurity density, E{r) varies by less than 50% over r [131 
Fig.l]. Thus it is appropriate to use the approximation 



Ve/h{r) = ^'^j;,E{r) , (2) 

where fif^f^ are the effective mobilities of the charge carriers. 

The differential equations, Ve,h(t) = drefi(t)/dt, for the radial drift, Te^hit), of 
the electrons and holes along a crystallographic axis can be solved. In the case 
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considered here, the electron drift dominates and any T dependence, fil^(T), 
is directly reflected in a T dependence of the rise time, t^-lT), which becomes 



fiT) ' 



(3) 



where C^et is a constant that was calculated from the detector geometry and 
its electrical parameters after solving the equations. For the detector under 
consideration it is Cdet = 5.9 ■ 10~^cm^/V. 

It is generally assumed that fJ^f^{T) can be described as 



where /xj is a constant, / is some functional form and T is given in Kelvin. 
Models based on the scattering of electrons and holes off the phonons in the 
lattice [H] suggest a form of /(T) = T~^/^. This is qualitatively supported by 
some measurements at higher impurity levels [T5| p. 108]. Other measurements, 
also at higher impurity levels, show a behavior of /(T) = T~^'^ or /(T) = 

The values of the electron mobility used in the simulation for the (110) axis 
at 77K can be translated into a value for /ig^(77K) of 1.4 ■ 10^cm^/(Vs). For 
/(T) = (T/K)-3/2^ the constant /if becomes /if = 9.4 ■ lO'^cm^ /(Vs). This 
provides an absolute model prediction to be compared to data. As it is not 
expected that the shapes of the pulses for core and segment vary signiflcantly 
within the T range, it also provides a prediction for 



5 Event selection 



Events induced by 122 keV photons were selected using the energies measured 
in the individual segments. Single segment events were selected by requiring 
energy deposits of less than 20 keV in all other segments. Figure H] shows 
the spectrum around 122 keV of the segment containing the (110) axis for 
single segment events from DS3 at Tmon=104K. The width of the peak was 
determined by a flt of a Gauss function and events within ±2a were selected. 
Typical signal to background ratios were 3:1 to 5:1 for DS2 and DS3 and 1:1 
to 2:1 for DSl. 



(4) 
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Energy [keV] 

Fig. 4. Energy spectrum around the 122 keV peak as measured by the segment on the 
(110) axis. Also shown is a fit with a Gauss function plus a first-order polynomial. 
The 2a limits are shown as dashed lines. 

6 Rise time determination 



The rise time of a measured pulse was not extracted directly because of the 
noise level observed during the measurements. Instead, a fitting procedure 
was used, where a simulated pulse was fitted to the measured pulse as a 
reference. This method was chosen, because it uses the full information of the 
measurement and is robust against noise which averages out. 




Charge [a.u.] 

Fig. 5. Typical noise distribution in arbitrary units. Also shown is a Gaussian fit. 
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Before the fit, tlie baseline was subtracted from each pulse. It was calculated 
as the average of the values recorded for the first 15 bins, 200 ns, of the pulse 
record. This was possible, because the DAQ was adjusted to record pulse 
information starting 200 ns before the trigger time. 

The noise level for each 15-minute interval was determined using the first 
200 ns of several hundred baseline subtracted pulses. Figure [5] shows a typical 
noise distribution. The unit of the charge is arbitrary. The noise is clearly 
Gaussian. Its level was typically about 3 % of the pulse amplitude. Each point 
of a measured pulse was assigned an uncertainty corresponding to the noise 
level of the dataset. 

The simulated pulses, C^m, were fitted separately to the measured core and 
segment pulses: 

^meas(^) = '^'^''^ ' C*sim (^/'^scale + '^offset)) (5) 

where the fitted parameters are the relative amplitudes, A'^'^, the time offsets 
of the measured pulses, tofjset) ^^e time scaling factors tgcaie- The important 
parameters are tgcaie which were converted into measured values of 

Only good fits were considered. A fit was qualified as good, if x^/ndf < xls- 
The values for xl s were chosen to reject 50 % of the pulses corresponding to 
the worst background levels observed. They were adjusted for each dataset 
and where 1.1(2.0) for DSl, 1.2(1.7) for DS2 and 1.7(1.9) for DS3 for the core 
(segment). Figures IH] and [7] show examples of two events with good and bad 
fits for both core and segment pulses. 

A bad fit could indicate that the interaction point was not close to the outer 
mantle. This can happen when a 122 keV photon penetrates deeper than in 
average or, if the event results from an interaction associated to the Compton 
background. The background also contains multiple interactions. 

The distributions of the fit parameters t^^'^jg for all good fits were fitted with 
Gauss functions and the resulting means, ^scaie-mean; converted to measured 
^c,s:io-90 values for a given position and temperature. 

Figure [8] shows the f^^aie distribution for the core pulses along the (110) axis 
at 77.4 K. The resulting tgcaie-mean ^^^Y close to 1 and results in a measured 
pulse length of = 287±2 ns. The uncertainty is statistical only. 

The measured value for the segment is t^''^^~^^ = 305± 2 ns. The agreement 
between simulation and measurement is excellent for the core while it is quite 
poor for the segments. The discrepancy might be explained by a combination 
of the poor understanding of the hole mobility and net charge carrier density 
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Fig. 6. Pulses in core ctnci segment of ci typicctl event from DS3 cit rmon — 

100 K. 

The fits with the simulated pulses are also shown. The fits are classified as good 
with x^/ndf = 1.06(1.05) for the core(segment). 

variations close to the core. The hole mobility determines the time when the 
10 % level is reached. The 90 % level could be influenced by a changing drift 
velocity of the electrons close to the core due to changing impurity levels. The 
segment pulses would be affected more than the core pulses due to the different 
strengths of weighting fields. The relative discrepancy should, however, not 
change within the temperature range considered. 



7 T dependence 

The values of t^,s:io-90 determined according to the procedure described in 
the previous section were used to study the temperature dependence of the 
rise time. The analysis was performed separately for the two crystallographic 
axes (110) and (100) using the datasets Al and A2 from Tabled! The drift 
along (100) is known to be faster at 77 K [9j. The T dependence is a priori 
assumed to be the same. 
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Fig. 7. Pulses in core and segment of a typical event from DS3 at Tmon = 100 K. 
The fits with the simulated pulses are also shown. The fits are classified as bad with 
X^/ndf = 1.87(2.29) for the core (segment). 




Fig. 8. Distribution of the scaling parameter f^^^^^ for measured core pulses alon^ 
the (110) at 77.4 K. Also shown is a Gaussian fit. 
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As explained in section [2l T is only known up to a measurement dependent 
shift AT for DS2 and DS3. When a functional form /, as introduced in Equa- 
tion m was fitted to rise-time data including DS2 or DS3, free parameters 
AT^''^ were introduced to the fit, respectively. 

Figure [9] shows the development of tf^^~^^ for the (110) axis together with the 
model prediction from section HI 

t,{T) = ^T'/' , (6) 
K 

and a fit with a simple Boltzmann-like ansatz 

tr{T)=po+pie-P'/^ , (7) 
where po, pi, p2 are free parameters. 

The model prediction does not describe the data at all. It correctly predicts 
the value of t^: 10-90 77.4 K, but it is clear that the predicted T dependence is 
not observed in the data. The decription cannot be improved by modifying the 
exponent. The more general ansatz U{T) = po + piT^^ can also not describe 
the data. The description is bad in the T range of DS2 and DS3 and completely 
misses the DSl point. However, the Boltzmann-like ansatz describes the data 
quite well. This is also true for the segment pulses as shown in Fig. [TDl For 
the segment pulses, the model prediction from section H] neither describes the 
T dependence nor the point at 77 K. The hole drift is more important in this 
case and as mentioned in section [3] the parameters for the hole drift are less 
well measured. 

The situation for the (100) axis is shown in Figs. [11] and [121 No data were 
available at 77 K. However, it is clear that a power law cannot describe the 
data for either core or segment, while a Boltzmann-like ansatz can be fitted 
quite well. 

The success of the Boltzmann-like ansatz could be connected to a change 
in the conductivity of the germanium crystal [17] in the T range considered. 
This might have affected the electric field strongly, because the detector was 
operated relatively closely to the full depletion voltage. 

Table [2] lists the values obtained in x^-fits for the parameters pi and p2 from 
Equation [7] which are strongly correlated. The values for p2 are also shown 
as energies, E, with p2 = E/{2k), where k is the Boltzmann constant. The 
uncertainties given are statistical only. The x^/ndf is typically around 5. 

The core and segment parameters agree within statistical uncertainties for 
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Fig. 9. Dependence of tf'^^^^^ on T along the (110) axis. Also shown are the model 
prediction described in section^] and a fit with a Boltzmann-like ansatz. Statistical 
uncertainties are smaller than the symbol size. 




Temperature [K] 

Fig. 10. Dependence oit^'^^^^^ on T along the (110) axis. Also shown are the model 
prediction described in sectionlHand a fit with a Boltzmann-like ansatz. Statistical 
uncertainties are shown, but are mostly smaller than the symbol size. 
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Fig. 11. Dependence of f^'^^^^^ on T along the (100) axis. Also shown are the model 
prediction described in section^] and a fit with a Boltzmann-like ansatz. Statistical 
uncertainties are smaller than the symbol size. 




Fig. 12. Dependence of t^'^^^^^ on T along the (100) axis. Also shown are the model 
prediction described in sectionlHand a fit with a Boltzmann-like ansatz. Statistical 
uncertainties are shown, but are mostly smaller than the symbol size. 
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Table 2 

Fit parameters determined for a Boltzmann-like ansatz using the datasets Al, A2, 
CI and C2 





orientation 


c, s 


Pi [-10^] 


P2 


E [meV] 


Al 




c 


17.7 ± 1.3 


913 ± 10 


158 ±2 


Al 






15.0 ± 2.3 


900 ± 20 


155 ±4 


A2 


(100) 


c 


3.0 ± 0.2 


630 ± 8 


109 ±2 


A2 


(100) 


s 


3.1 ± 0.4 


666 ± 18 


115±3 


CI 


(110) +5° 


c 


3.5 ± 0.3 


649 ± 12 


112±2 


CI 


(110) +5° 


s 


5.6 ± 0.7 


725 ± 18 


125 ±3 


C2 


(110) -15° 


c 


5.5 ± 0.5 


705 ± 13 


122 ±2 


C2 


(110) -15° 


s 


5.5 ± 0.9 


732 ± 22 


126 ±4 



the datasets Al and A2, respectively. However, the parameters for the two 
axes seem to be significantly different. A cross-check for the (110) axis, using 
dataset CI from Table [H did, however, not confirm this. The parameters 
should be the same, but they agree better with the parameters of the fit to 
A2. The inclusion of DS2 is what pushes the values of and p2 up in the fit to 
Al. As the absolute temperatures for DS2 and DS3 are not known, but have 
to be fitted, we used this cross-check to estimate the systematic uncertainty 
on E to be ?«40meV. The systematic uncertainty could also be connected 
to the misalignment present for DS3, even though such a strong deviation 
is not expected for a 5°offset. Another cross-check was performed using the 
dataset C2. At that position, the drift is not simply radial any more due to 
the distance from the crystallographic axes [ISl Fig.3], but the data could 
still be described by the Boltzmann-like ansatz. The fit results for C2 and CI 
are compatible. 

The fits were further analysed using the program package BAT p!8]. The val- 
ues listed in Table [2] were confirmed by BAT as the global mode. A detailed 
study of the probability densities revealed that the inclusion of DS2 causes 
elongated two-dimensional probablity densities with local minima, also point- 
ing to systematic uncertainties of up to !^40meV for the fits to Al and A2. 
The probability densities for the fits to CI and C2 are much more regular, 
pointing at systematic uncertainties of less than 20 meV. 

The drift along the two crystallographic axes were expected to have the same 
T dependence. The data, independent of fitting procedures, seem to indicate 
that the T dependence of the "slow axis", (110), is less pronounced. At some 
T, the rise time along the "fast axis", (100), seems to become as large as 
along the "slow axis". A definite conclusion would need better data. More 
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measurements at more angles and with an improved temperature monitoring 
are planned. 

It should be noted that while the rise times changed with temperature, the 
amplitudes remained stable, even to the highest temperature of about 130 K. 
This means that the detector could be operated stably at temperatures well 
above 100 K. 



8 Summary 



The temperature dependence of the rise time of pulses was studied between 
77 K and 130 K. The data cannot be described by the expected T~^/^ depen- 
dence of the charge carrier mobilities. Both core and segment pulses display a 
Boltzmann-like temperature dependence of the rise time. 

The significant temperature dependence of the rise time should be taken into 
account when pulses are simulated for specific experimental setups or used 
in pulse shape analyses. It should be noted that the detector used for the 
investigation could be stably operated up to temperatures of 130 K without 
an observable change in the energy scale. 
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